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ABSTRACT 

We present a sequence of toy models for irradiated planet atmospheres, in which the effects of ge- 
ometry and energy redistribution are modelled self-consistently. We use separate but coupled grey 
atmosphere models to treat the ingoing stellar irradiation and outgoing planetary reradiation. We 
investigate how observed quantities such as full phase secondary eclipses and orbital phase curves 
depend on various important parameters, such as the depth at which irradiation is absorbed and the 
depth at which energy is redistributed. We also compare our results to the more detailed radiative 
transfer models in the literature, in order to understand how those map onto the toy model param- 
eter space. Such an approach can prove complementary to more detailed calculations, in that they 
demonstrate, in a simple way, how the solutions change depending on where, and how, energy redistri- 
bution occurs. As an example of the value of such models, we demonstrate how energy redistribution 
and temperature equilibration at moderate optical depths can lead to temperature inversions in the 
planetary atmosphere, which may be of some relevance to recent observational findings. 

Subject headings: line: formation - radiative transfer ~ atmospheric effects - eclipses - planetary systems 



1. INTRODUCTION 

The past several years have seen a rapid and exciting increase in the amount of information available concerning 
the physical properties of extrasolar planets. In particular, the class of short orbital period planets, known informally 
as 'Hot Jupiters', has begun to yield information through a variety of observational channels. The discovery of 
systems in which the orbital plane is almost edge-on to the line of sight has resulted in the detection of planetary 
transits (Charbonneau et al. 2000; Henry et al. 2000; Bouchy et al. 2005). Optical and ultraviolet photometry 
and spectroscopy of such systems can yield information about the absorbative properties of the planet atmosphere 
• (Brown et al. 2001; Charbonneau et al. 2002; Vidal-Madjar et al. 2003; Vidal-Madjar et al. 2004). In addition, 
the unparalleled capabilities of the Spitzer Space Telescope has allowed the detection of thermal infrared emission 
(Charbonneau et al. 2005; Deming et al. 2005, Deming et al. 2006; Grillmair et al. 2007; Richardson et al. 2007) 
during secondary eclipses as well as the measurement of phase curves for non-transitting systems (Harrington et al. 
; 2006). 

This flurry of observational activity has spurred a similar level of activity on the theoretical front, with several 
independent models being used to infer the physical properties of the planets from the observations (Seager & Sasselov 
2000; Barman, Hauschildt & AUard 2001; Cho et al. 2003; Menou et al 2003; Sudarsky, Burrows & Hubeny 2003; 
Burrows, Hubeny & Sudarsky 2005; Dyudina et al. 2005; Cooper & Showman 2005; Seager et al. 2005; Barman, 
Hauschildt & AUard 2005; Fortney et al. 2005; Fortney 2005; Burrows, Sudarsky & Hubeny 2006; Cooper & Showman 
\ 2006; Fortney et al. 2006; Langton & Laughlin 2007). The resulting inferences are not always entirely in agreement. 
As an example, consider each group's response to the first detection of secondary eclipses. Burrows et al. (2005) 
concluded that the redistribution of energy from the substellar to the antistellar face of the planet was weak. On the 
other hand. Barman, Hauschildt & AUard (2005), as well as Fortney et al. (2005), concluded that the redistribution 
was strong, while the opinion of Seager et al. (2005) fell somewhere in between. Similarly, predictions for the effects 
of atmospheric fluid motions on the emergent intensity pattern differ substantially from group to group (Cho et al. 
2003; Cooper & Showman 2006; Langton & Laughlin 2007). 

This diversity of opinion is hardly surprising, given the complexity of the problem. The study of Hot Jupiter atmo- 
spheres requires that we understand the atmospheric chemistry, radiative transfer and hydrodynamics of the planetary 
atmosphere. These are affected by the thermal structure and evolution of the planet, by the effects of photochemistry 
induced by ultraviolet light from the star, the possible existence and behaviour of atmospheric condensates, whether 
they form clouds, and whether (or if) energy is redistributed across the surface of the planet or reradiated in situ. 
Furthermore, there are a variety of different observational probes, each of which provides information but in ways that 
differ, sometimes subtly, from other methods. The complexity of the coupled radiative transfer, chemical equilibrium 
and atmospheric flow models also restrict somewhat the Uexibility of the models in providing the true qualitative or 
geometric information necessary to develop a real understanding of what the observations are telling us. 

The overall intent in the pages to follow is pedagogical. Clearly, detailed model atmospheres are necessary to fully 
interpret the emerging observations. However, simplified models that focus on key physical elements can be useful to 
try and develop a qualitative understanding of what elements of the model each of the various extant observations is 
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actually probing. To do this, we construct a sequence of analytic and semi-analytic toy models. These models are 
chosen to be simple enough to be transparent while still capturing the qualitative features of the more detailed models. 
The hope is that one can use this approach to tailor our interpretations of the more detailed calculations to reach a 
proper synthesis between theory and the different observations that we have at our disposal. In § [2] we review our 
simplest model, based on a variation of the grey atmosphere model. In this first case we will assume no redistribution 
of energy across the face of the planet. In § [3] we then examine how the various available observations relate to the 
underlying structure of the model. In §|3]we then introduce the physics of energy redistribution into the model and 
examine how this affects the observations in § [5j Finally, in § [6l we place our models within the context of both 
observations and more detailed models in the literature, and examine possible implications for future generations of 
models. 

2. MODELS WITHOUT REDISTRIBUTION 

To begin, we briefly review the basic steps in calculating the grey atmosphere model. These may be found in almost 
any textbook on radiative transfer (e.g. Mihalas 1978) but it is worth reviewing them again as we will need to employ 
some subtle variations in subsequent sections. In doing so we will also derive a simple model which gives an explicit 
analytic description of the current philosophy in the literature for treating this situation. This will serve as a useful 
baseline for comparison in later sections. 

2.1. The Simplest Model 

In the grey approximation, we work with a frequency-integrated version of the radiative transfer equation 

l^^=I-S (1) 

where / is the specific intensity, and S is the source function. / is a function of optical depth t and angle cosine 
(relative to the normal direction) ^. In local thermodynamic equilibrium, where the optical depth is due to absorption 
rather than scattering, S = B, the frequency-integrated Planck function. Integration over angle yields the equation 
for flux conservation, which is satisfied as long as the conditions for radiative equilibrium are met, i.e. J = B. J is 
the mean intensity, the zeroth moment of / over angle. The first moment of equation ([T]) yields 

where H and K are the first and second moments of /. At this point, we adopt the Eddington approximation, assuming 
that J = 3K. This results in the equation 

3^^ = ^-5 = (3) 

which yields a solution for J(t), J = Jq + J\t. 

The constant Ji is determined from equation ([2]) and the lower boundary condition on the flux Fint = 4iJ entering the 
atmosphere from below. In the terminology of Mihalas (1978), F is the "Astrophysical Flux" and H is the "Eddington 
Flux" . Since S = B = J , one can also determine / from the formal solution to the first order differential equation ([T]) 



/(O, fi) = / [Jo + Jit] e-'l^- = Jo + -F,^t\x. (4) 

The final act is to determine the constant Jq by the calculation of the emergent flux at the surface (t — 0) from the 
expression for / and the condition that it must equal Fint 

/I ^1 ^0 

/^/(O, M)rfAi = 2 / iil{Q,ii)dpi + 2 \ pil{^,^i)dii. (5) 
-1 Jo J-i 

For most astronomical applications, the second term here is set to zero as it represents the contribution from rays 
directed inward at the surface, and most stars are not subject to any significant level of irradiation. However, in the 
case of Hot Jupiters, this particular term is of great importance as it introduces the irradiation from the parent star 
into the problem. The simple approximation used most often in the literature is uniform irradiation i.e. J(0, [x) = Jo 
for ij. < 0, although the reality of the situation is that all the incoming rays are coming from a single direction, that 
of the star. Nevertheless, if we proceed with the uniform irradiation approximation, we can use this to determine Jq 
and derive the final solution for the specific intensity 

J(O,M) = ^F„,i(^/. + 0+/Jo, (6) 

where we have included a dilution factor / in the spirit of the more sophisticated modellers. This is how the detailed 
models attempt to account for the difference between the true situation and the isotropic irradiation assumption. 
Because of the computation expense involved, these models are still calculated in the one dimensional case, but the 
incoming radiation is 'diluted' by / in order to approximately mimic the effect of averaging over the true non-uniform 
irradiation pattern. 
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An important point to note is that the net flux through any given point in the model atmosphere is stiU only the 
internally generated i^mt, since no energy from the irradiation is permanently absorbed by the planet in this baseline 
model - it is all reradiated. However, the flux the observer sees is the sum total of all flux travelling outwards from 
the planetary surface, and so includes the contribution of reradiated stellar light. In the case of a secondary eclipse, 
this is radiated from the surface in the direction of the star. (For this simple model the distinction is moot, but it will 
become important in subsequent sections). 

Ffuu =2 [ ^il{0, fi)dti = F,nt + flo- (7) 
Jo 

The traditional assumption found in the literature has been to adopt a value of / = 1/2 to represent the case of no 
energy redistribution over the surface of the planet. In this case the dilution is meant to represent the fact that the 
incoming Hux is absorbed from plane parallel rays spread over over a surface nRp (where Rp is the planet radius) but 
is reradiated isotropically from each point on the irradiated hemisphere and so the energy is emitted into a solid angle 
of 27r. If the energy is redistributed uniformly over the entire planet by fluid motions/winds and only then reradiated, 
then the angle is Att and the dilution factor is / = 1/4. 

2.2. A Better No- Redistribution Model 

The previous section gave an explicit example for how irradiated planets are currently treated using a one-dimensional 
model and a dilution factor. However, we can adapt this formalism easily to treat the geometry of the no-redistribution 
model properly. 

The first improvement to make is to dispense with the assumption of an isotropic irradiation. Consider instead the 
irradiation is due to a mono-directional beam^, so that /(0,/Lt) = IoS{9 + 9o)S{4> — (jio) for n < 0, and 9 and (j> are the 
poloidal and toroidal co-ordinates for a given point on the planetary surface. Furthermore, we can now choose & 
(jjQ SO that, for each point on the planet's substellar hemisphere, the irradiating beam is coming from the direction of 
the star. To achieve this, we simply note that /iq is the angle cosine between the present location and the substellar 
point. To determine the constant Jq in this case, we need to match to the outgoing flux Fint + A*o-^o- With this change, 
equation Q is replaced by 

/(O, fio,fJ^) = ^Fint (^^ + + t^olo (8) 

We omit / in this expression because it is no longer needed. This is now a solution that varies across the surface of 
the irradiated planet, as determined by the local angle of incidence of the irradiation. To obtain the flux we actually 
observe, we simply integrate over 9 and (j). Evaluating the integral at full phase (secondary eclipse), we find 

/■I 2 
Ffuii = 2 / /io/(0,/Zo,Aio)d/xo = Fi^t + -lo- (9) 

JO 

It is of interest to note that, in order to make equation ([7]) match this expression, we require / = |, rather than the 

value f = \ that is normally quoted. This is because our new model does not have a uniform distribution of flux over 
the surface of the planet - rather, it is hotter near the substellar point because that is where the irradiation is the 
most intense. Thus, the uniform dilution approximation underestimates the secondary eclipse flux. 

We can easily confirm global energy conservation in this model. With the monodirectional beam of irradiation, the 
total flux impinging on the surface of the planet is 

Firr^2 dfif) i ^iIoS{fi + ^o)dn] = Iq. (10) 







Using equation ([8|) to calculate the total emitted flux returns the same value. In the case of the isotropic irradiation 
model, the total flux emitted from the planet may be expressed as 



dfio [2 J fil{0,ti)dtij = 2F,nt + 2 flo (11) 

if we use equation ([6]). The equality of this with Firr is how we arrive at the value / = 1/2. 

Thus, our very first conclusion, even from a relatively simple model, is that the detailed models in the literature often 
dilute even their 'no-redistribution' models too much. Models for secondary eclipses should not be too democratically 
averaged over the surface of the planet, but should be weighted in favour of the hotter substellar point. This is 
in agreement with a similar comparison performed by Barman, Hauschildt & AUard (2005), who compared two 
dimensional no redistribution models with the same one dimensional models for / = 1/2 (their "a" is equivalent to 
our /). 

^ The readers who a re co ncerned about the use of the Eddington approximation in concert with a monodirectional beam should defer 
their skepticism until § 12.31 The discussion here is illustrative and this point will be addressed further when we come to the actual model 
we will use. 
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2.3. A Two Opacity Model 

We wish to improve our toy model still further, because there is an important piece of physics still missing. The 
wavelengths of the incoming stellar photons are very different from those of the reradiated photons. This is a con- 
sequence of the very different temperatures and hence spectral energy distributions of the star and the planet. The 
result is that the opacities and hence optical depths can be very different for the energy flow in the optical and infrared 
wavelength regimes. Thus, we will modify our model such that we treat the absorption of stellar irradiation seperately 
and with a different optical depth than we do the subsequent reradiation. 

In order to do this, let us treat the grey version of the radiative transfer equation in a 'two-stream' fashion, with 
two specific intensities, Ii to describe the optical radiation, and I2 to describe the infrared radiation. This allows us 
to treat the two streams with different opacities and emissivities. In this limit, we have two equations of radiative 
transfer, corresponding to the optical and infra-red respectively. 

f^^=S2-l2 (13) 

OT2 

where we have set Si, the source function in the optical, to zero, since we assume the planet has negligible emissivity in 
the optical. Furthermore, the constraint of radiative equilibrium requires that all the energy removed from the beam 
(in both wavelength regimes) is replaced by the infrared source function S2 , so that 

KiJi + K2J2 = K2S2, (14) 

where Ki and Ji are the appropriate opacities and mean intensities^. Thus our radiative equilibrium expression is 
slightly modified from the usual form, 

82 = J2 + —Jl= J2+lJl- (15) 

Within this approximation, the first moment equations become 

dHi 



dH2 



-Ji (16) 
^2 - J2. (17) 



(18) 



Using equation (jlSp. and ri — 7x2, we can combine these equations into 

dH2 _ dHi 

9t2 dT2 ' 

This amounts to the rather obvious statement that all the flux removed from the incoming beam in the optical must 
emerge in the infrared beam. 
The attenuation of the incident beam is described by 

/(T,M) = /(0,M)e^"/^ (19) 

where we have replaced T2 with r, our default optical depth scale henceforth. Recall that /i < 0, since this is an 
inwardly directed beam. The local flux (travelling inwards) is 

F= f 6{^-M f filoSifi + /io)e^"/''dM - -Moloe"^"/'^'' . (20) 



This is the flux flowing inwards due to the irradiation. We have thus a measure of how much energy is deposited 
at each depth. We assume the energy is absorbed, thermalised and now reradiated in the infrared. This is done by 
incorporating equation (I^Dl) into equation (|18p . 

^ = - J/oe--/-. (21) 

The gradient in H is negative, because the total energy flux carried outwards by the radiation fleld is less as one gets 
deeper into the planet. 

We once again make the Eddington approximation. This is a poor one for a monodirectional beam, but we note 
that the approximation is only applied here to the reradiation field. Once the energy is thermalised and included in 

^ This approach bears a superficial resemblance to the formalism of Hubeny, Burrows & Sudarsky (2003), in which the two opacities are 
derived from different frequency averages. However, that study depends critically on the frequency and temperature dependance of the 
opacity microphysics, as it was focussed on the possibility that multiple solutions might arise for a given physical situation. The multiple 
solutions derived by Hubeny et al. do not arise in our case, because we assume the opacity ratio 7 is a constant. 
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the reradiation energy flow, which is much more isotropically distributed, this is a more reasonable approximation. 
Thus we obtain the integrable equation 

= --nhe-^^'^" (22) 



which can be integrated to yield 



9t2 



4 7 



(23) 



The integration constant Ji is determined, as before, by the internal flux emerging from the planet (Ji = |-Fi„t). 
To determine the other constant, we need to ensure the total flux emerging at the top of the atmosphere reflects the 
appropriate sum of internal and absorbed flux. To that end, we calculate the specific intensity as a function of angle 
at the surface 



-f(0,Mo,M) 



dt 



Ji 







10 M 4-""r^ 4 7(1 + 

One may then integrate this over angle to determine the total flux emerging from the surface 



F 



^/(0,^)d/^ = Jo 



1 



3 /^o^o , 3 ^IqIq 



ln(l + ^) 



(24) 



(25) 



/o 2 2 72 2 7- fio 

This should match the sum of the internal radiation and flux incident on the top of the atmosphere from the star. 
Pint + Mo^Oi so that the final expression for the emergent specific intensity is 



7(0,^0, m) 



7 



(26) 



This expression now will allow us to calculate the emergent radiation field in a manner appropriate to any particular 
observation geometry, given a particular choice of observer orientation and orbital phase. The angle /io measures 
the angle for each point on the surface relative to the illuminating star, while fi measures the angle relative to the 
observer. The parameter Fint is the internal heat flux emerging from the planet itself, while Iq encodes the strength of 
the irradiation. The ratio of Iq relative to Fint will determine whether the observational appearance is dominated by 
the reradiation of absorbed energy. The one other factor in this equation is 7. This is the ratio of the optical depths 
for absorbed and reradiated energy. It can be tuned to reflect an optical depth for incoming stellar light that is either 
more (7 > 1) or less (7 < 1) opaque than that for the emerging thermal radiation from the planet. Note that this 
does not reduce to equation ([5]) in the case of 7 = 1, because the quantity / now refers only to infra-red light, and so 
the flux carried in the beam is not constant with depth (as it is in §[2?2])- 
In the deep absorption limit (7 — > 0), 



/(0,/io,/i) 



-F 



8 /io 



1 - 2/? 



0(7^ 



(27) 



which can be cast to zeroth order in terms of just an extra contribution to the internal flux, as one would expect 



/(0,/io,/i) 



(F.nt + Mo/o) M + X + ^7^0 (1 - 2^2) + 0(72). 



In the opposite (7 



r 

4 ' V" ■ 37 ■ 8 

00) limit (absorption high in the atmosphere), we have 



-Ft 



2 \l 



7 



Mo 



0(1/7^) 



(28) 



(29) 



The last term, despite the logarithmic function of 7, still falls off faster that 1/7^, so retaining only the leading term 
in 7 leaves 

2 

Mo ' 



7(0, Mo, m) 



Fint M + o + MO^O 



1 



7 



1 - 



1 



2^ 



0(ln7/7^). 



(30) 



So, for high altitude absorption, the zeroth order limit is the simple no-redistribution model of equation [H]* 

Thus, equation (|26p presents us with a simple, yet flexible model that we can use to try and understand how different 
potential observations probe the properties of Hot Jupiter atmospheres. 

3. SOME APPLICATIONS 

One of the primary goals of this study is to provide a flexible framework in which to understand how the observations 
probe different elements of the physics. So, let us now consider now how this model is reflected in the various kind of 
observations that are becoming available. 



Note that the apparent divergence of the second order term as ^ is misleading. One cannot maintain the hmit 7 - 
time as ^ ^ 0. One can verify this by taking the appropriate edge-on limit of the full expression, which remains finite. 



00 at the same 
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3.1. Secondary Eclipses 

A secondary eclipse occurs when a planet in an edge-on orbit (as measured by the observer) passes behind the 
star. The size of the flux decrement observed yields the flux emerging from the substellar side of the planet. More 
specifically, it represents the flux observed at full phase, when = /xq. Evaluating the emergent flux integral at this 
phase, we get 



1 



1 



1 

47 



1 



72 47(1 + 7) 



7 
2 



1 



1, 



■ 7 



Tn(l+7)-7-' ln(l + -) 



7 



(31) 

In principle, comparing the observed Ffuii to model values can tell us whether energy is being redistributed to the 
antistellar side, because then the observed value would be less than expected from the no-redistribution model. We 
can express this as an effective value of / by dividing the reradiation term by /q, so that 



feff - o 



1 1 

1 + ^ + ^- 



47 72 47(1 -f 7) 



J-7'- AM1 + 7)-7'M1 + -) 



7 



(32) 



Despite the algebraic length of the expression, there is actually remarkably little variation as a function of 7. As 
7 — > 0, the limit is 



17 3 

h -7 - 

24 8 



0(7^) 



while the limit as 7 ^ cxd is 



/■ 



^ + ^ + 0(ln7/7=^) 



(33) 



(34) 



Figure [T] shows the function feffil)- There is only a change of about 6% between high and low values of 7. We may 
thus conclude that, in case of no energy redistribution, it matters little where in the atmosphere the energy is actually 
absorbed (at least for secondary eclipse measurements). We find that / = 2/3 in the large 7 limit, and in the low 7 
limit, / = 17/24 = 0.708. We note that this latter value corresponds to the same full-phase fraction derived by Milne 
(1926) in his treatment of the 'reflection effect' in close binaries. His treatment amounts to a special case (7 = 1) of 
the above treatment. The insensitivity of / to 7 is encouraging, because it suggests that the comparison of models 
to secondary eclipse fluxes is robust, in the sense that any flux deficit is likely to indeed be due to the process of 
redistribution, and not strongly affected by whether the energy is absorbed high in the atmosphere or deep. 

3.2. Limb Darkening 

Of course, there is much more information potentially available than just the flux at a single orientation. To start, 
let us consider the planetary limb-darkening law. Neglecting the internal flux term, at full phase we have 

-1 3 MO \ 3 { tio\'^ ^ { t^o^^ ^■n(^ 



-^(0,Mo.Mo) 
7(0,1,1) 



Mo- 



-) 



1 



4 7(1+7) + 27 



3 1 



|:Fln(l+7) 



(35) 



which is shown in Figure[2]for three values of 7 = 0.1,1 &10. Also shown is the normal limb-darkening law for internally 
generated heat 

^(0,Mo) 3 / 2\ 

7(o:Ty = 5r^3j' 

and the limb-darkening law derived by Milne (1926) in the two-beam approximation 

/(O, ^0, Aio) _ / 1 + 2/io 



7(0,1,1) 



3 



(37) 



We see that all the irradiation-driven profiles are much more limb-darkened than the traditional internally generated 
profile. This is to be expected, since, near the limb, the irradiation is being absorbed at oblique angles and so spread 
over more area. Absorption high in the atmosphere makes the limb even darker, but the effect is small relative to the 
dominant effect of the irradiation geometry. The steep dependance on /iq is another way to illustrate why the true 
value of / is / = 2/3 rather than / = 1/2; the region near the substellar point, which is also the hottest, contributes 
disproportionately to the reradiated flux observed at full phase. 



3.3. Secondary Transit Shapes 

During secondary transit, the face of the planet is steadily eclipsed by the star. Thus, the manner in which the 
total light from the system decreases and increases can, in principle, tell us how light is distributed on the face of the 
planet. Figure [3] shows the shape of the secondary transit ingress (i.e. the drop in the system flux as the star occults 
the planet) for three cases. In each case, we assume the occulting body is a sharp, vertical boundary. The solid line 
indicates a planet with the no-redistribution illumination pattern oriented exactly edge-on. The dotted line is for the 
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same model but now with an inclination of 2 degrees. The dashed line indicates a uniform hemisphere planet oriented 
exactly edge-on. Each case assumed the planet radius was 0.1 times that of the star and that the orbital semi-major 
axis is 12 stellar radii. Separating out inclination effects is relatively easy based on the duration of the eclipse, but 
the measuring differences in illumination pattern by determining the shape of the ingress is going to be quite difficult. 
The differences are only a few percent at most. This is consistent with other studies concerned with using the shape 
of secondary transit to infer the presence of non-uniform brightness (Williams et al. 2006; Rauscher et al. 2007a). 
Greater potential for constraining models is likely to be found in detecting variations in the timing (Williams et al. 
2006) or amplitude (Rauscher et al. 2007b) of maximum depth. 

Overall then, it appears as though the properties of secondary eclipses are not strongly affected by the value of 
7 in this model. This is encouraging, since that means we can make robust inferences about the nature of energy 
redistribution. 

3.4. Phase Curves 

Another way to probe the asymmetry of the flux distribution emerging from the planet is to measure the phase curve 
- the variation in the observed brightness of the planet as it orbits the star. For the simplest case, let us consider the 
phase curve for a planet which has two hemispheres, each of uniform brightness but at different levels. This amounts 
to a model in the spirit of the default model in § 12.11 Furthermore, let us set the antistellar side to zero brightness, 
since it is the difference in flux between sub- and antistellar sides that is important. 

A location on the surface of the sphere is identified by a unit vector defined in terms of two angles, 9 and </), 

n — (cos (j) sin 6, sin cj) sin 9, cos 9) . (38) 

We consider the direction towards the observer to lie along the x axis, so that the angle cosine between the observer 
and local normal on the surface of the sphere is then fi = cos sin 0. We also need to specify the direction from which 
the irradiation comes, using the unit vector 

n, = (cos A cos i, sin A, — cos A sin i) . (39) 

where the angle A can circulate in a plane rotated by an angle i about the y-axis. Thus, the angle cosine between this 
direction and local normal is 

/io = n.n» = cos (j) cos A sin 9 cos i + sin (j) sin A cos 9 — cos 9 cos A sin i (40) 

To take account of the phase, we note that we need not account for the true orbital motion, but only the rotation of the 
irradiation pattern, since the actual orbital displacement is unresolved by any applicable photometric measurement. 
Thus, our phase curve is described by the circulation of A. For the edge-on case {i = 0) this has an analytic solution 
for the phase curve 

/A+7r/2 pTT 
d(t> / sin^ 6* COS </>d6' (41) 
-tt/2 Jo 

for A < 0. This applies for brightness of unity for fiQ = cos{(j) — A) sm9 > and zero otherwise. Evaluation of the 
integral yields 

^ = |[l + cosA] (42) 

which also holds for A > 0. Figure S] shows the phase curves for this model and a range of inclinations. 

We can repeat the same exercise for the model in § 12. 2[ except that the surface brightness is no longer uniform, since 
there is an extra factor of /iq in the integrand in equation (I41|) . The solution is then 

[A cos A - sin A] (A < 0) 

[A cos A -sin A] (A > 0). (43) 

The lower panel of Figure [4] shows the appropriate phase curves. The qualitative difference is that this model has 
broader troughs and narrower peaks than the uniform hemisphere model. This is a consequence of the concentration 
of the illumination near the substellar point, so that the phase curve is more sensitive to how this particular region 
is projected onto the line of sight. Some readers will no doubt note that this is the same as the phase curve for 
a Lambert sphere (e.g. Russell 1916). Although the Lambert sphere is usually invoked to describe reflected light, 
the conditions under which it is derived (single, isotropic scattering) amount to the same conditions as applicable to 
thermal reradiation. 

In Figure [5] we show a comparison between these and the phase curves of Barman et al (2005), calculated with a 
full radiative transfer model for the no-redistribution case. The models show excellent agreement. This suggests that, 
for reasonable models, the phase curve is determined primarily by the geometry and less by the effects of changing 
temperature and atmospheric chemistry towards the limb. Indeed, we do not show any phase curves indicating the 
effect of varying the 7 parameter, because the variation is quite small (^ 7% change in amplitude between 7 = 10^^ 
and 7 = 100, with very little change in shape). 
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3.5. Temperature Profile 

We can also calculate the atmospheric temperature structure of these models, using the mean intensity J and the 
condition of radiative equilibrium, J — aT'^ /tt. 

\- T A 

3 



1/4 



eff 



1 



1 + J^) _ 3mo^_^,/,.„ 



Mo 



4 7 



(44) 



where i*i„t = aT^ff. Ideally, we would like to frame this using a quantity more physically useful than Iq, so we shall 



define an equivalent temperature by Jo 



rpi 



eff 



Reframing the equation then, we have 
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Figure [6] and Figure [7] show the temperature profiles at a range of angles mo, for two cases (7 = 0.1 and 7 = 10), and 
assuming (To/Te//)* = 100. For the case where energy is absorbed high in the atmosphere the temperature profiles 
look pretty similar to the unirradiated case, although the asymptotic value of the temperature as r — > is higher the 
closer one gets to the substellar point. In the case where energy is mostly absorbed below the reradiative photosphere 
we find that the temperature profile has two regions where temperature is almost constant - the traditional one at 
low optical depth but also a second one at r > 1, where the bulk of the incoming flux is absorbed. This is similar 
in character to the profiles found in many detailed radiative transfer models (e.g., Seager, Whitney & Sasselov 2000; 
Barman et al 2005, Burrows, Sudarsky & Hubeny 2006). 

The origin of this intermediate plateau is easily understood. The linearity of the equations of radiative transfer 
and equilibrium means that the above solution corresponds to the superposition of the solutions to two independent 
problems ~ the standard solution radiating Fint at the boundary with no irradiation plus the solution corresponding to 
zero net energy transfer with the irradiation boundary condition imposed. As noted by Milne, the temperature profile 
arising from this latter solution consists of a gentle transition between two asymptotic temperatures (at high and low 
optical depth). The difference between the two asymptotes is driven by the last term in equation (|45|) . leading to a 
temperature jump AT/Tq — (3/ig/47)^/^, and so is larger for smaller 7. Thus, in Figure [7] we see the behaviour of the 
zero flux solution at low optical depth and then the transition to the normal thermal profile deeper inside the planet. 
The same behaviour is found in Figure [6] but the difference between the two asymptotic temperatures is smaller (larger 
7) and can barely be discerned in the figure. 

3.6. Spectral Lines 

The presence or absence of spectral lines in secondary eclipse measurements is a topic of recent excitement (Grillmair 
et al. 2007; Richardson et al. 2007; Swain et al. 2007). Although our model is founded on the grey atmosphere 
approximation, we can make simple models of spectral line formation by considering a simple line formation model 
imposed on the continuum structure represented by our model. 

Consider a spectral absorption line resulting from an excess opacity at some specific wavelength. The line opacity at 
this wavelength is represented by r' = (1 + a)T. We may thus determine the strength of the line using the absorption 
depth A=l- Fiine/F^onti where 
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S{T)E2{T)dT 

is the continuum flux, S is the standard source function, and 

/>C30 

Fune=2 Sl,ne{T')E2{T')dT' 
JQ 

is the flux at line center. The function E2{x) — t^^e~^*dt. 

We assume local LTE, so that S = B = T'^ , and we use the background temperature structure from equation 
The resulting integration yields the expected Fcont — T^ft + Mo^^o ■ The same calculation to determine Fu 
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that we replace t in equation 
the limit Te// ^ To, is 



45|l with r'/(l + a) before evaluating the integral. The resulting expression for A, in 
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Thus, we can expect strong absorption lines in the case where the optical opacity is much less than the opacity in 
the infrared, but weak lines when the optical opacity is large. Similar considerations apply in the case of emission 
lines (a < and A —A). Such behaviour is also consistent with the model proposed in Richardson et al. (2007) 
for example, which invoke TiO or Silicate clouds as sources which can absorb energy at low pressures (high in the 
atmosphere), in order to explain the lack of observable water absorption. 

3.7. Transmission Spectroscopy 

Spectral line formation is also relevant to the question of absorption line spectroscopy during primary transit (Char- 
bonneau et al. 2002; Vidal-Madjar et al. 2003, 2004). However, it is somewhat different in detail, since the observation 
is essentially measuring the attenuation of a parallel beam illuminating and (partially) passing through the limb of the 
planet. Furthermore, we need to make a further approximation, since the calculation of the transverse optical depth 
at the limb requires a relation between the physical length scale and the optical depth scale. 

We make use of the fact that the absorption occurs high in the atmosphere, where the temperature profile is 
essentially isothermal. We note, however, that the r — > asymptote is still a function of angle from the substellar 
point. This can be seen from Figure [S] or [71 Assuming that the equation of state is that of an ideal gas, we can 
approximate the pressure as P = pgH, where H{iio) is the local scale height. We can then use hydrostatic equilibrium 
to derive a relation between density and optical depth, so that p ~ t/{kH{pq)), where t is still the same local radial 
optical depth as before. We can then also infer the relationship between optical depth and radial distance r 

r = Toe"'-/^(''°\ (51) 

where tq is a normalisation constant than can be chosen to fit this approximation (only valid high in the atmosphere) 
to the true relation. Note that this relationship is a function of angle, so that a surface of constant optical depth is 
not spherical. 

Turning to the optical depth for absorption through the limb, which we denote as r^, we have 

R ff(Aio) (r2 - 

where R is the cylindrical radius and so + — r^ and z = fior (We also include a factor 7 to relate this optical 
absorption to the infrared opacity represented by k). We can make the evaluation of this integral easier by noting 
that the scale height is small compared to the radial and cylindrical scales for the atmospheric base. As such, the 
exponential in the integrand will limit significant contributions to scales such that r — R < Hq, where ~ H{fiQ = 0) 
is the atmospheric scale height at the limb. This comes about because the density drops rapidly as the light path 
moves away from /xq = and results in a negligible contribution to the absorption. With this approximation the 
integral becomes 

r, = jro e-^/"" y-'/'e-ydy = jto e-^/^^ . (53) 

We can see then, that at cylindrical radii R, the optical depth to transmission is 7(27ri?/iJo)^^^ larger than the radial 
optical depth at that location, meaning that the ti — 1 surface will lie higher in the atmosphere, where r < 1 (because 
i?> ifo). 

Using equation (|53p. we can estimate the so-called 'transit radius effect' (Baraffe et al. 2003, Burrows et al. 2003), 
which is the difference between the radius measured by observing a planetary transit, and that calculated in a traditional 
stellar model (which is usually determined at some fixed radial optical depth Tf). Setting = 1, to determine Rt, 
and setting t — Tf to determine Rt — AR, we infer 







T£ = 27 / Kpdz = 27T0 / — — ^—^ „ (52) 



Ai?=^ln 
2 



27ri?T 



(54) 



Thus, the transit radius increment is proportional to the local scale height at the limb, and only logarithmically 
dependant on 7 or r^. Recall that Hq oc T^ff/g, so that there is an inverse relationship with gravity. More surprisingly, 
there is no direct dependance on the intensity of the illumination, since the asymptotic temperature at the limb {fxo —>- 0) 
does not depend on Tq. Of course, the irradiation does slow the planet cooling, so that Te// would be higher for the 
same planet under strong illumination, but this is an indirect effect. 



4. ENERGY REDISTRIBUTION 

The models in the previous section considered the absorption and re-emission of all radiation locally, i.e. there 
was no energy transfer between different locations on the planet. However, we have abundant evidence for strong 
atmospheric flows in our own giant planets and so it is quite possible that winds and atmospheric circulation could 
redistribute the absorbed energy over a significant fraction of the planet before it is reradiated. This would have a 
significant influence on the observations, because the no redistribution model leads to large temperature differences 
between the substellar and antistellar sides of the planet. 
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4.1. Simplest model 

Thus, we turn now to the case of redistribution of energy over the surface of the planet. In this section, we will return 
to the approximation of § 12.21 wherein we impose the effects of irradiation by use of the outer boundary condition only. 
We do this in order to isolate the physics of redistribution separately from any considerations about where the energy 
is actually deposited. In the next section, we will combine our redistribution model with the more detailed 7 models 
outlined in § 12.31 Furthermore, we will assume the redistribution of energy over the surface happens at a single, fixed 
optical depth Tw This introduces a new source/sink term in the radiative transfer equation to reflect the fact that 
energy may be added or removed from the net flux, depending on location. Equation ([3]) is thereby replaced by 

1 52 J 

-g^ ^ -ASir -T^) (55) 

where A is a function (of position on the planet) describing the local energy input/loss due to atmospheric redistribution. 

The solution to this equation is of the standard form J — Jq + Jit except that now the first derivative experiences 
a discontinuous jump at r = r^. Taking into account the internal flux Fint we find 

J^Jo+ (^^F.r^t + 3A^ T (r < T^) (56) 
3 

= Jo + ^Fmt (r>T^) (57) 

As before, we determine the constant Jq by calculating first 7(0, /i) and using that to calculate the net flux emerging 
from the surface 

/I 
/i/(0,M)d/^ 2 / fiI{0,fi)dfi + 2 fil{0,fi)dii (58) 
-1 Jo J-i 

= -/io/o + 2 / til{0,fi)dn (59) 
Jo 

where we have again used as our upper boundary condition for incoming radiation a model in which the incoming 
radiation is modelled as a beam of intensity /q coming from a direction iiq. Furthermore, this net flux has to equal 
the sum of F^t and the effect of the A term, whether it is positive or negative. Thus, F = F^t + 4A = F^t + Ai^. 
Using this to determine Jo results in a final answer for the specific intensity of 



3 / 2\ 3 
7(0, Ai) = -F,nt ( Ai + 3 j +Aio/o + ^AF 



(60) 



where 7?„(x) = e "dt is the standard function found in radiative transfer theory and we have made use of the 
identity xEn-i = — {n — l)En. In the case of A = 0, this reverts to equation ([8]). 

To model the redistribution, we need to specify the form of the function A. We implement redistribution by imposing 
temperature equilibration on curves imposed by the fluid circulation. Recall that, in LTE, the specific intensity J — B, 
which is the local energy density, so that the redistribution of energy essentially involves equilibration of the mean 
intensity J around the appropriate curve. In particular, using equation ([56]) . we can calculate J^, the mean intensity 
at the optical depth where redistribution is assumed to occur 



3 / 2\ 3 
Jw = -^F^nt ( + 3 ) + l^olo + ^ AF 



^ + 2 (e-^™ - 27;4(t^)) 



(61) 



To complete the model, we also need to specify the geometry, i.e., along which curve is the temperature equilibrated? 
There are a variety of models for the fluid circulation on the surfaces of irradiated planets (see § 16. 3[) , and each offers 
a different possibility. We will base our model on the simulations of Cho et al (2003), who find that the strong winds 
maintain a banded structure, with a hot equator and cold poles. This is qualitatively similar to the handedness of 
Jupiter and other solar system planets, although the bands are fewer and wider because of the different rotational state 
of the system. Thus, we enforce equilibration of on bands of constant 9p (the angle measured from the pole of the 
planet). Recall that the angle cosine iiq reflects the angle between the local normal and the direction of the irradiating 
star. Casting this in terms of angles relative to the pole, yields /io = sin^pcos^p. The second term in equation (|6ip is 
linear in fiQ and so we wish A to have a similar dependance in order to cancel out some or all of that angular variation, 
in the interests of redistribution over 0p for constant 9p. Thus, we set A = Pi + P2 cos0p on the substellar side, and 
A = /3i on the antistellar side (where we do not need the (j)p dependance to cancel out any irradiation term). Energy 
conservation requires that, if we remove energy at one point on the curve, we must inject it again somewhere else, so 
that J^^ Ad(/)p = 0, and so 7r/3i + /32 = 0.. Furthermore total equilibration requires that all terms in equation (j6ip that 
depend on (j)p must cancel, yielding an expression for Pi 

/3i - T^^^ (62) 
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and the resulting mean intensity is 

The intensity at the surface then foUows, 
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(63) 



(64) 



(65) 



This quantity e encapsulates how the surface emission is affected by the depth at which the redistribution occurs. In the 
limit oo, e ^ 1, in which case equation ((64)) becomes equation dS]). This is because, if the redistribution occurs 

deep in the atmosphere, none of the irradiation actually reaches that depth, and the absorbed energy is reradiated in 
situ as before. On the other hand, as Tu, — > 0, e — > Tu,/2 — > 0, and the mean intensity distribution is independent of 
(jjp, although not 0p, reflecting our approximation that the temperature structure is now banded. 
Let us now examine what fraction of the total luminosity is reradiated on the substellar side and the antistellar side. 



Integrating F = 2 M^(0, ^o, IJ')dn, and then integrating over 9p and 0p, we calculate 
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(67) 



We see that Lanti as — > oo, as expected. As — > on the other hand, Lanti asymptotes to almost Lg^^^ 
as expected in the case of complete redistribution. We will postpone calculation of most of the observables until 
§ [SI but we will calculate the expected full-phase flux, as this is the easiest way to estimate whether redistribution is 
taking place. Recall that to calculate the full-phase flux, we calculate Fjuii = 2 yLto/(0, /zq, /xo)d/j,o, i.e. we impose 
the condition ii = fiQ. However, to do this easily, we need to change angular co-ordinate systems. The equilibration 
was performed on contours that encircle the pole of the planet, but the integration at full- phase is most conveniently 
performed in angles measured from the substellar point. Thus, let us define the angle 4>s as that which circulates about 
the substellar point, and 6s the angle relative to that point. As a result, sin6'p = (1 — sin^ 9s sin^ 0s) 1/2. The angle 
^0 = cos 6s and we can thus perform the integral over 6s and (ps . Figure [5] shows the resulting flux, expressed in terms 
of the effective /e// parameter of equation ([7]). We see that the value tends to the expected no-redistribution value of 
2/3 as —> oo and comes close to the expected full redistribution value of 0.25 as — > 0. The actual asymptotic 
value is 
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(68) 
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The resulting phase curves are shown in Figure O For deep redistribution, the phase curve is very similar to the no 
redistribution curve, while redistribution at low optical depth suggests almost complete uniformity. 

4.2. Redistribution in the j-model 

We now wish to combine the effects of our models in 12.31 and S 14.11 that is, allow for both a different optical depth 
in the optical and IR and for redistribution. This will constitute our final model. Thus, equation ([3]) now has two 
source/sink terms 

^ = -^lloe--'^^'-'' " 3A5(t - T^). (69) 
We follow the same approach as before to derive the mean intensity at 
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Once again, we wish to equilibrate this quantity on contours of constant 6p. The functional dependance on 0p is now 
more complicated than a single cosine, so let us denote the last term as Iq sin0p<I>(7, cos(/)p), so that all the variation 
is contained in the function $. Thus, equilibration requires 



AF = /3i + /32$ 



(71) 
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to cancel out the variations on the substeUar side, while preserving the integral over the contour. Thus 

2 ^ losmOp ( 1 /■'^/^ 



Pi = 7. 9 7 - / ^ d(t)r,\ (72) 

_ (73) 



37r + f + 2(e-^". -2S4(t^)) 



We will henceforth denote 



<$>-(-/ $ ■ (74) 



-7r/2 



The resulting specific intensity at the surface is 



, 3 / 2\ /osin6'c /osin6'„, , 3uo^o 



(75) 



7 1 + 7M/Mo. 

where e is defined as in equation (j64p . This looks much the same as the simpler version, with an extra term on the 
end that accounts for the fact that the energy is deposited over a range in optical depths, and some may be deposited 
above Tw and some below. 

We can now calculate how much of the reradiated luminosity is emitted from the substeUar side, as a function of 7 
and Figure [TOl shows the substeUar and antistellar fractions as a function of r^, for three different cases, 7 = 100, 
7 = 1 and 7 = 0.01. At large 7, we see the expected behaviour, where large r^, implies that very little energy is 
reradiated on the antistellar side and small Tw leads to almost complete equilibration. For 7 3> 1, $ — > cos (^^ and the 
last term in equation ((75l) is ~ 0(7"^) even if — > 0, so that we recover the same equations as in 12.11 This tells 
us that absorption high in the atmosphere can still result in complete equilibration as long as is relatively small. 
However, absorption high in the atmosphere results in little equilibration if r^, > 1. 

In the case for 7=1 we see qualitatively similar behaviour, except that the solution does not asymptote to complete 
equality as r^, — > 0. The difference is even more extreme for 7 <C 1. Here we see that there is an intermediate regime 
ijw 3) where almost complete equality is established, but then that solution diverges again as — > 0, with fully 
three quarters of the luminosity emerging on the substeUar side. This is because, in the limit 7 — > (corresponding to 
deep deposition of energy), the limit of equation ((75|) is 

/(0,^o,M,7 0, ^ 0) = + ( -^/osin6'pcos0p ) /x, (76) 



27r2 V47r 

in which the second term now behaves like a traditional flux contribution, powered by irradiation energy deposited well 
below the photosphere. The physical origin of this behaviour is that, in this limit, many photons do not contribute to 
the temperature at low r because they are absorbed deep in the atmosphere (because the optical opacity is so low) 
and then escape from below t„ (which is small in this limit). 

We can summarise the various effects by examining Figure 111! which shows contours indicating the fraction of the 
total absorbed luminosity that is reradiated on the substeUar side. We see that there are two regimes which allow 
for approximate equality between substeUar and antistellar sides. The first is for 7^1 and r^, <C 1. This regime 
corresponds to the case when both absorption and redistribution both occur high in the atmosphere, as might occur 
if one has a layer of stratospheric clouds and strong stratospheric winds. The second is the perhaps more traditional 
case of 7 <C 1 and 1 < r^, < 10, where energy deposition and energy redistribution occur at moderate (infra-red) 
optical depths. Note also that the contours tend to change direction near the line 7 ~ t^. Above this line much of 
the irradiation energy is absorbed in regions of the atmosphere above the redistribution region, while below it the 
incoming photons penetrate to regions below before being absorbed. 

5. SOME MORE APPLICATIONS 

Let us now revisit the various observational probes with our new models and examine how things are changed by 
the inclusion of redistribution. 

5.1. Secondary Eclipses 

We can, once again, calculate the full phase flux, as observed during secondary eclipse, using equation ([75]) Once 
again, the flux integral is performed with ^ = /iq. The result is shown in Figure [121 where we have again cast the result 
in terms of the /e// parameter. We see that, for ^ 1, we recover the no- redistribution limit for all 7 as expected. 
For a broad range of low and intermediate values of t^, the value of /e// varies between 0.4 and 0.5. For low 7, the 
value of feff actually has a minimum near Tw ~ 1. For large 7, the behaviour is more like the simple model, with near 
complete redistribution at low t^, changing monotonically to no redistribution at large r^,. 

The diversity of values obtained demonstrates the complexity of inferring planetary properties from as simple an 
observation as comparing a flux level at a single phase to the models. We note that the value /e// — 0.25 that one 
expects from the naive full redistribution model is realised only in the r^, ^ 1, 7 ;3> 1 limit and that /e// attains 
intermediate values ~ 0.4-0.6 over most of the parameter space. Thus, while we may be able to demonstrate, using 
this kind of observation, that a planet possesses some level of redistribution, it will be difficult to determine specific 
parameters without additional information. 
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5.2. Phase Curves 

This behaviour is also reflected in the diversity of phase curves, shown in Figure[T21 The three panels show the phase 
curves (for edge-on orbits) for three choices of 7. In each panel we show curves for three different r^,. As expected, 
large r^, makes for a strong variation with phase, while low r^o results in a very flat phase curve. At low 7 values, this 
distinction is blurred somewhat. 

Figure [14] summarises this behaviour. We show here the amplitude of the observable phase variation 5f, from peak 
to trough, expressed once again in terms of feff- We see, not surprisingly, similar behaviour to Figure In the 
large limit, we recover Sf > 0.6, essentially no redistribution. On the other hand, amplitudes of Sf < 0.1 are 
obtainable in the same two regimes as identified above. Once again, intermediate behaviour is recovered over most of 
the parameter space, with phase variations ~ 50% of the no-redistribution model expected. 

5.3. Limb Darkening 

Figure [15] shows the influence of redistribution on the limb-darkening of the planet, when observed at full phase. With 
the model of redistribution adopted here, the limb is no longer at a constant intensity, but rather shows a variation 
oc sin^p. Figure [TSl shows contours of the contrast observed at full phase between the intensity at 9p ~ tt/2, <j)p = tt/2 
(the limb at the equator) and the intensity at the center dp = tt/2, (j>p = 0. We see that deep redistribution (t„ ^ 1) 
yields low values, indicating strongly darkened limbs, as in the no-redistribution case. We also find dark limbs in 
the limit 7 <C 1 and <C 1, where the energy is absorbed deep in the star but redistribution occurs high in the 
atmosphere. However, for 7^1 (absorption high in the atmosphere), we actually find equality and possibly even 
limb-brightening. Similar effects have been seen in more detailed models, for instance when the addition of metals 
results in a large amount of stratospheric absorption (Fortney et al 2006) . 

5.4. Temperature Profiles 
We can also examine the temperature profile. For t < t^, 
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For T > Tu], the temperature profile still obeys equation (|l5|) . 

Figure [16] shows temperature profiles for two cases where 7 = 100. In the upper panel, we show the result for 
shallow redistribution (r^, — 0.1) and in the lower panel, deep redistribution (tw = 10). In each case, the solid curve 
is the temperature profile at the substellar point {9p — tt/2, (pp — 0) and the dotted curve is for the limb {9p = tt/2, 
(pp = tt/2). In each case, the ratio {Tq/T,,/ f)'^ = 100 and the temperature is measured relative to the true internal 
Teff. In the case of = 10, the effect of the temperature jumps are largely erased by the time r — > 0, presumably 
because most of the absorbed energy never reaches r = 10 when 7 = 100. In this case, a temperature inversion is 
evident, and an important part of the energy balance at moderate r is thus the transport of energy inwards by thermal 
radiation from overlying layers heated by the stellar photons. Thus, equilibration at r^, =0.1 does a much better job 
of smoothing out the temperature variations. 

In the opposite limit (7 = 0.01), deep absorption, the temperature profiles are shown in Figure [17] In this case, 
the opposite behaviour is seen. The deep redistribution case is effective, because most of the irradiation energy is 
deposited in deeper layers and is included in the thermal balance at t = t^,. In the shallow redistribution case, the 
temperatures near the substellar point are much hotter at r ~ 1 and so this case shows a strong asymmetry. There 
is an interesting center-to-limb temperature inversion at low r for this case, but it will be difficult to probe with any 
vertically integrated measures. 

5.5. Spectral Lines 

We can evaluate the integrals for spectral line formation in a similar fashion as in § 13.61 although the algebra is 
somewhat more painful this time, and so we will not write the full expression down here. Figure [TS] shows the contours 
of A = 1 — Fiine/Fcont that rcsult from the appropriate integration in the case where a = 0.3, for the usual range of 7 
and . The calculation is once again at full phase orientation. We see that the strength of the line is relatively weak 
for 7 > 1 and rather insensitive to . This is what we expect on the basis of Figure 1161 because the temperature 
profile is close to isothermal. For 7 < 1, the results are more sensitive to t^, and can be quite strong in the case 
where t„, has a value only slightly larger than unity. This is because the combination of sharp temperature jumps and 
different opacity scales can lead to large contrasts between line and continuum. 

6. DISCUSSION 

The overall goal of this exercise is to develop a simple framework within which one can examine the qualitative 
features of the irradiated atmosphere problem. Clearly, a simple, semi-analytic model such as this cannot replace 
numerical models that incorporate detailed treatments of radiative transfer, chemical equilibrium (and disequilibrium) 
and hydrodynamics. However, the complexity of the various physical inputs into these models often results in a 
situation in which the physical understanding lags behind the calculational output. As a result, simplified models 
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provide a pedagogical framework which can guide the more detailed models. This is likely to be particularly useful 
as we try to develop physical models for the redistribution of energy within an irradiated planet atmosphere. This is 
likely to be a difficult and costly enterprise, and the simplified models may help to illuminate the nature of the desired 
solution and possible pitfalls on the path to their realisation. As a first step, we wish to cast our models in terms of 
the more detailed models in the literature, in the sense of describing which of our simple models correspond to the 
more detailed models. 

The first attempts to treat the structure of irradiated giant planets (Saumon et al. 1996; Guillot et al. 1996) used 
essentially grey atmosphere models with a simple asymptotic temperature criterion. The first real attempt to model 
the radiative transfer under conditions appropriate to Hot Jupiters was by Seager & Sasselov (1998). These models 
were dust free, and the incoming optical/UV radiation was found to penetrate well below the infra-red photosphere. 
These models therefore fall into the 7 < 1 category, and establish a theme for most of the subsequent literature. This 
was followed by Barman, Hauschildt & AUard (2001), who produced two sets of models, "clear", in which condensates 
were removed by a rainout process, or "cloudy", in which the condensates were assumed to form clouds in situ. The 
cloudy models were found to yield clouds extending from the photosphere upwards to low optical depths, and most 
likely correspond to the 7 > 1 case. This is supported by the claim that this led to much smoother emergent spectra 
than in the clear CBtSG, clS expected for spectral features in the 7^1 model. 

Sudarsky, Burrows & Hubeny (2003) present a variety of models, covering a large range of irradiations. The models 
they classify as classes IV and V are most appropriate to Hot Jupiters. The majority of their models clearly fall within 
the 7 < 1 regime. This can best be illustrated by the temperature profiles. Their Figure 12 shows how the temperature 
profile of a planet changes as it is moved closer to a GOV primary. At distances < 1 AU, the temperature profile 
developes a plateau at pressures >1 bar, evolving to a shape very similar to that in Figure [HI For the hottest planets. 
Iron and Forsterite condensates are possible at high altitude, leading to smaller plateaus, as might be expected for 
larger 7 atmospheres. The clouds in this case are localised, however, and so the analogy with the 7 models is not as 
good as for the Barman et al. (2001) models. 

For the purposes of calibration the models of Barman, Hauschildt & AUard (2005) are of particular interest because 
they calculated a truly two-dimensional model and so produce temperature profiles as a function of angle from the 
substellar point, rather than just a single 'representative' profile. We can thus accurately calibrate the effective 7 of 
their model by making use of the temperature value of the observed plateau in the model. We calibrate Tq using the 
known properties of the HD 209458b system, Tq = T^f f{R*l af/"^ = 2030 K. Taking the profile from their Figure 4 
for HD 209458b at the substellar point, the high temperature plateau has a value ^ 2500 K. We identify this value 
with the T ^ 00 limit of the irradiation term in equation ([45]). We use this expression for clarity since the Barman 
et al. models are for the case of no redistribution, but equation ([77)1 reduces to this in the limit of r^u 00. The 
temperature of this plateau is 
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We can solve this to infer 7 ~ 0.45 for /io 1. Figure [T9l shows the expression 
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as a function of 7. We see that the inner plateau temperature increases as 7 decreases, with the energy getting 
deposited at greater depth. The lower asymptote for these values of 7 and Tq is derived from T — 3/iQ/47 (shown 
as the lower curve in Figure [T9|) . and yields a temperature ~ 1800 K. Furthermore, Barman et al. use an internal 
luminosity effective temperature of T^^jf — 230 K. The resulting temperature profiles are shown in Figure BUI The 
agreement with Figure 4 of Barman et al is excellent for r > 1 (corresponding to P > 10 mbar). At low optical depths, 
the temperature drops lower in the detailed models, although it does eventually asymptote to higher temperatures off 
the edge of the given plot (Barman, private communication). This behaviour presumably reflects the non-monotonic 
nature of the real opacity at high altitudes. A general feature of these kinds of model is that (infrared) optical depth 
unity lies somewhere in the middle of the transition between the two plateau. This can be seen in other models as 
well, such as from the location of the 'brightness temperature region' in Figure 1 of Fortney et al. (2005). 



6.1. Redistribution as dilution 

The detection of secondary eclipses in the HD209458 and TrES-1 systems generated a flurry of theoretical activity. Of 
particular interest was the question of whether energy was transferred from the day side to the night side. However, all 
of these models were flux-conserving and so approximate redistribution by diluting the incoming flux by an appropriate 
factor /. Burrows, Sudarsky & Hubeny (2005) compare their / = 0.25 (full redistribution) models to the observed 
data and find that they are too faint, suggesting that less dilution is best. BSH also express surprise that adding 
iron or forsterite clouds to their models does not change the predicted fluxes by an appreciable amount. We can 
understand this by noting that increasing 7 at fixed r^, ~ 100 does not change the predicted /e//, as was noted in 
§ 13.11 and can be seen from Figure [121 Fortney et al. (2005) find that they cannot fit the two IRAC measurements 
for TrES-1 with a model without redistribution, but are able to almost fit the data with a diluted radiation field 
appropriate to significant redistribution. The temperature profiles in their Figure 1 suggest again that their model lies 
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in the 7 < 1 regime, with the characteristic two-plateau behaviour very evident. They note that the inclusion of an 
arbitrary extra continuum opacity in the range 2-6/j,m would improve the fit, and which would also move their model 
to larger 7. Seager et al. (2005) include, as an additional constraint, a limit on the flux variation at wavelengths 
~ 2.2/im (Richardson et al. 2003). This places limits on the depth of observable water bands at near-IR wavelengths. 
Combining this with the secondary eclipse measurements, they found that the clear (7 < 1) models were ruled out, 
but models with thick, absorbative clouds did match the models. In this case clouds composed of Silicates extend 
from 0.8-5 mbar and Iron condensates from 1.3-10 mbar. From their Figure 4, we can infer directly the approximate 
7 values appropriate to their models. In the case of the cloud-free models, the optical photosphere lies at greater 
pressures than the infrared photosphere, and 7 0.1, while the cloudy model has 7 ~ 10. Similarly Barman et al. 
(2005) find that their no-redistribution model overpredicts the flux at IRAC wavelengths, so they too infer that some 
level of energy redistribution must be occuring. 

Subsequently, Fortney et al. (2006) have used enhanced metallicity models to predict the structure of the newly 
discovered hot Jupiter HD149026 b. Their models with 3 times solar metallicity, when calculated in the '47r' (full 
redistribution) case, have little cloud absorption because the relevant condensation curves lie at pressures of hundreds 
of bars. In the '27r' (no redistribution) case, the models result in hot stratospheres because the TiO and VO absorption 
is present at lower pressures for such hot models. The comparison between the resulting temperature profiles can also 
be understood in the context of our models, since the increase in 7 for the no redistribution model leads to a smaller 
temperature difference between the two plateau, as observed. In this case though, even the hotter model still lies in 
the 7 < 1 regime. 

Burrows, Sudarsky & Hubeny (2006) make a first attempt at predicting phase curves while including a primitive 
redistribution prescription, in which they remove a fraction P„ of the irradiation energy (i.e. this is still a dilution 
model) and then incorporate it as an extra internal energy source on the antistellar side. Their temperature profiles 
again very clearly conform to the 7 < 1 structure we have seen before, consistent with their neglect of clouds. They 
do, however, discuss one 'extreme' model with significant heating of the stratosphere using significant TiO and VO 
absorption in the upper atmosphere. This may indeed lie with the 7 > 1 regime. Burrows et al present phase curves 
for a variety of values of P„ and bandpasses, and note that even significant redistribution can still leave some remnant 
phase variation because of the non-uniform temperature pattern on the substellar side. This has some qualitative 
resemblance to the kinds of variation seen in Figure UM 

6.2. Including redistribution in the radiation field directly 

However, all of the previous attempts to model redistribution have relied on some kind of dilution of the original 
radiation, rather than removal of the energy at the appropriate point in the atmosphere. The models discussed here 
allow us to examine how the temperature profile changes if we actually include sources and sinks at different points in 
the atmosphere. 

Let us consider the effect of redistribution at different r^,, using our calibrated version of the Barman et al. (2005) 
model, shown in Figure [20l Figure [211 shows the effect, on the temperature profile at the substellar point, of moving 
to progressively smaller values. Because this is for the substellar point, redistribution removes energy at the specified 
location, resulting in some pronounced local temperature inversions. In fact, for ~ 10^-10^, there is a pronounced 
local peak in the temperature, resulting from a competition between the local deposition of irradiation energy and 
the removal by redistribution. At low enough r^, we start to see a significant drop in the stratospheric temperature. 
The equatorial variation in temperature profile, from center to limb, is shown in Figure [22l for the case r^, = 100. We 
see that the effect of redistribution changes from a sink to a source for fiQ < 0.4, although the general character of 
a hot substellar point is retained, albeit with reduced amplitude. In the case where r^, = 1, shown in Figure [23l the 
temperature at high altitudes is much more uniform. 

These effects apply for the case of 7 < 1, which applies to most of the models in the literature. However, there 
is tentative evidence that some systems may favour the 7 > 1 regime, so let us consider redistribution in this case. 
First, we show in Figure [Ml the effect of changing 7, while keeping all other inputs fixed. The dotted line shows our 
calibrated model for 7 — 0.45. Lower values of 7 result in even higher temperature plateau, as expected, while larger 
7 drives the atmosphere closer to an isothermal state. This is precisely why the observed lack of spectral lines favours 
a high 7 model. Figure [25] then shows the effects of different t^^ at the substellar point, but this time for the 7 = 10 
model. The temperature jumps arc less in this case, but there are still temperature inversions directly above the layer 
where redistribution occurs. However, the stratospheric temperatures do return to close to their value in the Tw 00 
case, with large differences only becoming apparent at t^, < 5. We note, however, that it is precisely this region (large 
7, moderate r^,) that seems to be preferred by the data for HD209458b and HD189733b. 

6.3. Hydrodynamic redistribution models 

We have thus far discussed the attempts to model energy redistribution in terms of calculations of radiative transfer. 
A second, and mostly independent, body of work is focussed on the hydrodynamics of irradiated atmospheres and how 
this might, in turn, influence the redistribution of energy. 

The complexity of the physics involved has resulted in somewhat divergent treatments of the phenomenon, each 
placing emphasis on different aspects of the phenomenon. Showman & Guillot (2002) performed three-dimensional fluid 
dynamical simulations of the atmosphere, flnding the existence of broad east-west jets and strong residual temperature 
contrasts at the photosphere. Their models also suggested that some fraction of the irradiation energy might get 
transported deeper into the planet, where it might help to explain the anomalously large radius of HD 209458b. Iro 
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et al. (2005), however, found that very httle energy was deposited at pressures > 100 bars. Iro et al. use a coarse 
version of a one-dimensional radiative transfer model, and include temporally variable forcing to mimic the planetary 
rotation. However, they do confirm the existence of a strong phase variation within the context of their model, which 
also fits into the 7 < 1 class of our models. In an updated version of their three-dimensional model. Cooper & 
Showman (2005) quantify their phase variation, predicting that the equatorial fluid motions should shift the peak of 
the phase variation ~ 60° downstream, although subsequent revisions of the model quoted in Knutson et al. (2007b) 
find a smaller phase shift. Similar behaviour is found by Langton & Laughlin (2007). An entirely different approach 
was taken by Cho et al. (2003; 2008), using a high resolution, two-dimensional model of the fluid equations on a 
rotating sphere, treated in the shallow water approximation. Cho et al. find generic behaviour in which the principal 
hydrodynamical feature is a small number of broad equatorial bands, accompanied by coherent polar vortices that 
rotate about the pole on timescales of several planetary rotation periods. The motion of the temperature field is found 
to be quite diverse, depending on the strength of the heating and momentum transport in the atmosphere and where 
the energy is deposited. 

Fluid motions are responsible not only for equilibration, but potentially also mixing. Cooper & Showman (2006) 
showed that hydrodynamic mixing processes can smooth out temperature and compositional differences between 
different parts of the atmosphere. They predict that this will lead to homogeneous Carbon chemistry, with most of 
the Carbon residing in CO and very little in CII4 at low (<1 bar) pressures. Furthermore, the resulting temperature 
structure is considerably more isothermal and model spectra show weak features (Fortney et al. 2006), in similar 
fashion to the observations. It is interesting to consider these latter models in the light of our 7 parameter. A 
uniform vertical redistribution of energy by hydrodynamic mixing can be considered equivalent to uniform absorption 
throughout the atmosphere, and can lead to a higher effective 7 than might be inferred simply from the opacity. 

Examination of Fortney et al. (2006) suggests another application of our method. They use a global dynamical 
simulation of the atmosphere to calculate how the temperature profile varies across the face of the planet and then 
calculate many one-dimensional radiative transfer models to quantify the likely spectral appearance as determined by 
the true two-dimensional surface variations. However, they were not able to iterate between the dynamical model and 
radiative transfer model to achieve a true radiative equilibrium because of the prohibitive computational expense that 
would entail. A model such as the one described here, suitably parameterised by an effective 7, could be used as an 
intermediate approximation to enforce a global energetic equilibrium while still retaining the spirit of the radiative 
transfer model. The resulting temperature distribution could then be used to calculate the detailed appearance. 

7. CONCLUSION 

In summary, the intent of the models presented here was to understand the qualitative behaviour of radiative transfer 
in hot Jupiter atmospheres. The simplified nature of the models is dictated by the difficulties of including both full 
radiative transfer and detailed hydrodynamics, both of which arc important parts of a proper understanding of the 
phenomenon. In particular we have studied how the models change depending on the parameters 7 (which refiects 
the difference in opacity for incoming optical-wavelength radiation and outgoing infrared-wavelength radiation) and 
Tw (the depth at which energy is redistributed horizontally across the planets atmosphere). 

We have also studied how the various new observational probes of Hot Jupiter atmospheres are influenced by 
changes in these parameters. We find that models in which redistribution occurs at large infra-red optical depths 
redistribute little energy to the night side (as expected). Perhaps more surprisingly, we find that significant day/night 
flux differences can also be found in the case Tw < 1, as long as 7 < 1 as well. We find that day-side temperature 
inversions are a generic feature of atmospheres in which the redistribution occurs at moderate (r^ ~ 1-1000) optical 
depths. This is of interest because of recent claims that temperature inversions have been found in HD 209458b 
(Knutson et al. 2007). The initial explanation for such behaviour is high altitude absorption (Burrows et al. 2007; 
Fortney et al. 2007), although the very low albedo (Rowe et al. 2007) restricts the absorbing species to be a poor 
reflector. Some level of redistribution is required to explain the HD 209458 b secondary eclipse amplitude and (lack of) 
phase variations, so that another potential explanation for the temperature inversion may simply be that it is related 
to the depth at which the energy is removed from the day-side. Certainly, figures [22l and [23l show that some parts of 
the planet experience pronounced temperature inversions, even without a significant high altitude absorber. However, 
it will require detailed radiative transfer models to confirm whether this is a viable model or not. 

Our initial exploration has been far from exhaustive. One important advantage of these toy grey atmosphere models 
is their fiexibility, and we have studied how energy redistribution changes the character of the models. This will 
hopefully provide a guide as to how one might perform similar calculations using the more detailed models, which 
presently can only consider changes at either the top or bottom of the atmosphere. Furthermore, we have also not 
exhausted the parameter space of toy models. For instance, one can imagine many other redistribution schemes apart 
from the one we have implemented here. 

This work has been supported in part by NASA ATP contract NNG04GK53G and by funds associated with Spitzer 
Space Telescope program GO-20101. The author thanks Sara Seager and Travis Barman for helpful comments. 
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Fig. 1. — The flux observed at superior conjunction, as a function of the energy deposition depth parameter 7, and expressed as an 
effective value of the parameter /. We see that the valeue of /e/ / is rather insensitive to 7. 
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Fig. 2. — The dotted line indicates the normal limb darkening law for an unirradiated grey atmosphere. The three solid lines indicate 
the same for our irradiated no-redistribution models, for three different values of the deposition depth parameter 7. The values are (from 
the left) 7 = 10, 7 = 1 and 7 = 0.1. The dashed line represents the expression in Milne (1926), equation 14. 
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Fig. 3. — The solid and dotted curves represent the secondary ecUpse ingress for the no-redistribution model, for zero and 2° inclination 
respectively. The dashed line is the edge-on ingress for a uniform hemisphere model. 




Fig. 4. — The upper panel shows phase curves for the case of the uniformly bright hemisphere model. The lower panel shows phase 
curves for the case of the no-redistribution mode l. C urves are shown in each case for inclinations of 0° ,45° ,60° and 75° . The edge-on cases 
are given by the analytic solutions in equations I I42I I and II43II respectively. 
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Fig. 5. — The solid and open circles represent phase curves for 24fj,m and Sfim respectively, calculate d us ing the models of Barman, 
Hauschildt & AUard (2005). The dotted line represents the no-redistribution model described in equation II43I I. The solid line is the same 
solution, but superimposed on a DC level of 0.19 on this scale, to represent a contribution from internal luminosity. The excellent agreement 
suggests that the phase curves are determined primarily by geometrical effects. 




Fig. 6. — The temperature profiles for angles of 0,30,60 and 90° to the line of sight, for the case 7 = 0.1, i.e. where the irradiation 
penetrates to below the reradiation photosphere. We have also used (To/Teyy)'' = 100. 




Fig. 7. — The temperature profiles for angles of 0,30,60 and 90° to the line of sight, for the case 7 = 10, i.e. where the irradiation is 
absorbed high in the atmosphere. We have also used (To/T^jj)* = 100. 
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Fig. 8. — The solid curve shows the effective / value for the full phase flux in our complete redistribution model, as a function of the 
redistribution depth Tu,. The dotted lines indicate the canonical values of 0.25, 0.5 and 2/3 discussed in the text. 
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Fig. 9. — The solid line represents the phase curve for a planet in an edge-on orbit in whieh the redistribution occurs at Tyj = 1. The 
dotted and dashed lines indicate the same, but for Tw = 0.1 and t„; = 10. One can observe the transition from almost no redistribution to 
almost complete uniformity as r^, decreases. The flux level is parameterised in terms of effective /. 
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Fig. 10. — The three panels show how the relative luminosity of the substellar and antistellar sides vary with Tw, for three diflferent 
values of 7. For large 7 (upper panel), we see a smooth transition from complete redistribution at low Tu, to essentially no redistribution 
at large r^,. For lower values of 7 however, the dependance is more eomplicated. 




Fig. 11. — The solid lines are contours of isub/^totaii ranging from 0.6 to 0.9. The dashed contour is a value of 0.55. These latter 
contours enclose that region of the phase spa<;e where the energy is approximately evenly redistributed. 




Fig. 12. — The solid contours indicate values of the full-phase scahng parameter / that correspond to canonical values 1/2 and 2/3. 
Dashed curves indicate other values as labelled. 
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Fig. 13. — In each panel, corresponding to a different value of 7 as labelled, the solid curve indicates Tw = 100, the dotted curve Tu, = 1 
and the dashed curve Tui = 0.01. In the middle panel, the dotted and dashed curves lie very close to one another. The nature of the phase 
curve can be seen to depend clearly on both 7 and Tyj. 
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Fig. 14. — The contours plotted here arc for the quantity Sf, which is the difference in feff when observed at full phase and zero phase. 
We have assumed an edge-on orbit here. The dotted and dashed curves represent 5f = 0.6 and Sf = 0.1 respectively. The solid contours 
lie between these two values, in increments of 0.1. We see that phase variations of moderateamplitude are expected over a wide range of 
parameters. 
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Fig. 15. — The contours indicate ratios of the intensity at the equator hmb to the intensity at the substellar point, when the planet is 
observed at full phase. The solid contours indicate values of 0.5 and 1.0, while the dashed contours arc 0.25 and 0.75. We sec that deep 
redistribution indicates strong limb darkening, while high altitude absorption with moderately shallow redistribution can actually lead to 
limb brightening. 
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Fig. 16. — The solid lines show the temperature profile for the substellar point in each case. The dotted lines indicate the temperature 
profile at the limb (for 6p = it/2). The upper panel indicates redistribution is efficient in the case for low T and large 7 but inefficient for 
large t and large 7. 
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Fig. 17. — The solid lines show the temperature profile for the substellar point in each case. The dotted lines indicate the temperature 
profile at the limb (for 6p = 7r/2). The upper panel indicates redistribution is inefficient in the case for low r and low 7 but efficient for 
large r and small 7. 




Fig. 18. — The solid contours indicate values of spectral line strength A. We see that the strongest lines occur for 7 <C 1 and Tw ~ 2. 
For 7 > 1, the lines are relatively weak and show little dependance on Tw 
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Fig. 19. — The upper solid line is the function T(7). For 7 <g 1, the function asymptotes to — > 1 + 3/27, indicated by the dotted line. 
The dashed line corresponds to the solution for the Barman et al. (2005) case. The lower solid line is the corresponding function for the 
lower temperature plateau in our simple model. 
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Fig. 20. — The solid curves show the temperature profiles for = 1.0,0.8,0.6,0.4,0.2 and (from top to bottom). The value of 7 = 0.45 
was chosen to match the normalisation of tlic models of Barman ct al. (2005). The other parameters To and T^ff were chosen based on 
the values used by Barman et al. The horizontal axis is really r, but we have scaled this down by a factor 63 to mimic pressure in bars. 
(We expect the scaling to be linear for a plane parallel atmosphere with constant opacity. In fact, using the known gravity of HD 209458b, 
we can even derive the effective mean IR opacity ~ 0.07cm? . 
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Fig. 21. — The curves show how the temperature profile of the 7 = 0.45 model at the substellar point changes as we decrease t^, the 
optical depth of the redistribution layer. Each curve is labelled with the appropriate value of Tw The effect of redistribution at the hottest 
point on the star is clearly to remove energy and transport it to cooler parts of the star. Thus, we see depressions in the temperature 
profile near r^, . However, for values of Tw ~ lO'^-lO*, we see that this region lies in the middle of the temperature plateau where much of 
the energy is being absorbed. The combination of these two effects leads to conspicuous temperature bumps. 
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Fig. 22. — The curves show how the temperature profile of the 7 = 0.45 model varies along the equator from eeiiter to limb, assuming 
T-w = 100. The curves are for values of /iq running from 1.0 (top) to 0.0 (bottom), in increments of 0.2. The heavy solid line is the profile 
at the limb, which also applies on the nonirradiated antistellar side of the planet. We see that the removal of energy by irradiation near 
the center changes into an injection of energy near the limb and on the antistellaj: side. 
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Fig. 23. — The curves show how the temperature profile of the 7 = 0.45 model varies along the equator from eenter to limb, assuming 
= 1. The curves are for values of fiQ running from 1.0 (top) to 0.0 (bottom), in increments of 0.2. The heavy solid line is the profile 
at the limb, which also applies on the nonirradiated antistellar side of the planet. We see that the removal of energy by irradiation near 
the center changes into an injection of energy near the limb and on the antistellaj: side, with a net result of enforcing a very uniform 
temperature at lower pressures. 
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Fig. 24. — The curves show the temperature profiles at the substellar point {fio = 1) as we increase 7 from a value of 0.01 to 7 = 10. 
As 7 increases, the incoming energy is absorbed at lower and lower pressures, and the temper atur e jump decreases, eventually leading to 
an almost isothermal profile. The dotted line is for the value 7 = 0.45, as calibrated for Figure [20l 
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Fig. 25. — Redistribution in the 7 = 10 model is not quite as dramatic, although we still find temperature inversions, which get deeper 
as Tw decreases. Recall that 7 ~ 10 and r^, ~ 3 is a good model for the HD209458b case, suggesting that this atmosphere may have a 
temperature inversion at the 20 millibar level. 



